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Summary 


Global asymptotic stability of a class of nonlinear multibody flexible space-structures 
under dissipative compensation is established. Two cases are considered. The first case 
allows unlimited nonlinear motion of the entire system and uses quaternion feedback. 

The second case assumes that the central body motion is in the linear range although the 
other bodies can undergo unrestricted nonlinear motion. For both cases, the stability is 
proved to be robust to inherent nonlinearities and modeling uncertainties. Furthermore, 
for the second case, the stability is also shown to be robust to certain actuator and sensor 
nonlinearities. The stability proofs use the Lyapunov approach and exploit the inherent 
passivity of such systems. The results are applicable to a wide class of systems, including 
flexible space-structures with articulated flexible appendages. 

Introduction 

Many space missions envisioned for the future will require flexible multibody space 
systems such as space platforms with multiple articulated payloads, and space-based 
manipulators used for satellite assembly and servicing. Such systems are expected to 
have significant flexibility in the structural members as well as joints. Control systems 
design for such systems is a difficult problem because of the highly nonlinear dynamics, 
large number of significant elastic modes with low inherent damping, and uncertainties in 
the mathematical model. The literature discussed below contains a number of important 
stability results for certain subclasses of this problem; e.g., linear flexible structures, 
nonlinear multibody rigid structures, and most recently, multibody flexible structures. 
Under certain conditions, the input-output maps for such systems can be shown to be 
“passive” [1]. The Lyapunov and passivity approaches are used in [2] to demonstrate global 
asymptotic stability of linear flexible space structures (with no articulated appendages) 
for a class of dissipative compensators. The stability properties were shown to be robust 
to first-order actuator dynamics and certain actuator/sensor nonlinearities. Multibody 
rigid structures comprise another class of systems for which stability results have been 
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advanced. Ideally, subject to certain restrictions, these systems can be categorized as 
“natural systems” [3]. Such systems are known to exhibit global asymptotic stability under 
proportional-and-derivative (PD) control. Upon recognition that rigid manipulators belong 
to the class of natural systems, a number of researchers [4], [5], [6], [7] etc., have established 
global asymptotic stability of terrestrial rigid manipulators employing PD control with 
gravity compensation. Stability of tracking controllers was investigated in [8] and [9] for 
rigid manipulators. In [10], an extension of the results of [9] to the exponentially stable 
tracking control for flexible multilink manipulators, local to the desired trajectory, was 
obtained. Lyapunov stability of multilink flexible systems was addressed in [11], However, 
the global asymptotic stability for nonlinear, multilink, flexible space-structures has not 
been addressed in the literature, and that is the subject of this paper. 

We consider a complete nonlinear rotational dynamic model of a multibody flexible 
spacecraft which is assumed to have a branched geometry, i.e., it has a central flexible body 
to which various flexible appendage bodies are attached (Figure 1). Global asymptotic 
stability of such systems controlled by dissipative controllers is proved. In many applications 
the central body has a large mass and moments of inertia as compared to any other 
appendage bodies. For this case, the effects of realistic nonlinearities in the actuators and 
sensors are investigated when the central body is in the attitude-hold configuration. The 
proofs given use the Lyapunov approach. For systems with linear actuators and sensors, 
the stability proof by Lyapunov’s method can take a simpler form if the Work-Energy 
Rate principle [11] is used. However, since the Work-Energy Rate principle is applicable 
only when the system is holonomic and scleronomic in nature, we have used a more direct 
approach in evaluating the time derivative of the Lyapunov function so that the results are 
more general. 

Symbols 

B control influence matrix 

C Coriolis and centrifugal force matrix 
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G p , G p 

G r 
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M(p) 
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P, P 
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X, X 

vp 

Vr 
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a 
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d t 
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damping matrix 

structural damping matrix 

position gain matrices 

rate gain matrix 

stiffness matrix of the system 

stiffness matrix for the flexible degrees of freedom 

number of rigid body degrees of freedom 

Lagrangian of the system 

mass-inertia matrix of the system 

number of total degrees of freedom 

generalized coordinate vectors 

vector of rigid body coordinates 

vector of flexural coordinates 

skew-symmetric matrix 

vector of control input 

Lyapunov function candidate 

state vectors 

position output 

rate output 

state vector 

the quaternion 

ith component of quaternion 

vector part of quaternion 

ith component of unit vector along eigen axis 

scalar defined as (a 4 - 1) 

integral of w 

Euler angle vector 
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ip a i actuator nonlinearity (ith loop) 

rp pt position sensor nonlinearity (ith loop) 

Vvi rate sensor nonlinearity (ith loop) 

u angular velocity vector for the central body 

u>x skew-symmetric matrix of the components of u; 

0 vector of rotational degrees of freedom between rigid bodies 

Mathematical Model 

The class of systems considered here have a branched configuration a 5 shown in Figure 1. 
Each branch by itself could be a serial multibody structure. For the sake of simplicity, and 
without loss of generality, we will consider a spacecraft with only one such branch (Figure 2) 
where each appendage body has one degree of freedom (hinge) with respect to the previous 
body in the chain. The results obtained in this paper, however, will also be applicable to 
the general case with multiple branches. Consider the spacecraft consisting of a central 
flexible body and a chain of ( k - 3) flexible links. The central body has three rigid rotational 
degrees of freedom, and each link is connected by one rotational degree of freedom to the 
neighboring link. That means there are k rigid body degrees of freedom. The Lagrangian for 
the system under consideration can be given by 

L-p T M(p)p-q T Kq ( 1 ) 

where, p = {u T ,0 T ,q T } T ; w is the 3 x 1 inertial angular velocity vector for the central body; 

0 = ( 0 \, 02 , where 6 t denotes the joint angle for the ith joint expressed in body- 

fixed coordinates; q is the (n - k) vector of flexible degrees of freedom (modal amplitudes); 
M(p) = M T (p) > 0 is the configuration-dependent mass-inertia matrix, and K is the 
symmetric positive definite stiffness matrix related to the flexible degrees of freedom. Using 
the Lagrangian (1) the following equations of motion are obtained. The details of the 
derivation of math model can be found in [13]. 

M(p)p + C{p,p)p+ Dp + I<p - B t u (2) 
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where {p} = {7 T ,0 T ,q T } T , and 7 = w. C(p,p) corresponds to Coriolis and centrifugal forces; D 
is the symmetric, positive semidefinite damping matrix; B = [I kx k 0 kx(n-k) ] is the control 
influence matrix and u is the ^-vector of applied torques. The first three components of u 
are the torques applied to the central body by attitude control actuators (one about each 
body-fixed axis), and the remaining components are the torques applied at the (k - 3) joints. 
K and D are symmetric, positive semidefinite stiffness and damping matrices: 


K = 


Ofcxfc 

.^(n-k)xk 


^kx(n-k) 
K (n-k)x[n-k) . 


D = 


0 kxk 

.°(n-fc)xfc 


Ofcx(n-fc) 
D(n—k)x(n — k) . 


(3) 


where K and D are symmetric positive definite. The angular measurements for the central 
body are Euler angles (not the vector 7), whereas remaining angular measurements between 
bodies are relative angles. One important inherent property (which we shall call “property 
5”) of such systems that is crucial to the stability results to be presented is given below. 
Property 5 : For the system represented by equation (2), the matrix (5 M - C) is skew- 
symmetric. 

Outline of Proof : Using the indicial notation, the it, j-th element of C(p,p ) is defined as 


1 dM k j dM k , dM t] . 

c « * = E + ~ "sr lp ' 

t=l 1=1 J 

Similarly, the kj- th component of the time derivative of the inertia matrix, M(p), is given by 
the chain rule as 


"*i- = E-5T» 


i=1 


Now if we define the matrix S = ( \M - C), then the Jfc, j-th element of S is given by 


S k] = i^Mkj ~ C kj ) 

_ 1 A \ dMkj_ , dM kj dM kt _ dMjj ' 

~ 2 h\ i dpi 1 dp, dpj dp k \ 

dMij 6M ki 1 


P, 


= £ 


Pi 


d Pk dpj \ 

Since the inertia matrix is symmetric, i.e., Mjj = My, it follows from above equation by 
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interchanging the indices k and j that 


Sjk = -Skj 

Which means the matrix 5 is skew-symmetric. ■ 

It is assumed that the sensors consist of angular position and rate sensors which are 
collocated with the torque actuators. The sensor outputs are then given by: 

y p = Bp and y r = Bp (4) 

where p = (jj t , 9 t , q T ) T wherein 77 is the Euler angle vector for the central body. y P = (v T ,o T ) T 
and y T = (u> T 1 0 T ) T are measured angular position and rate vectors, respectively. It is 
assumed that the body rate measurements, cj, are available via rate gyros. 

Quaternion as a Measure of Attitude 

The orientation of a free-floating body can be minimally represented by a 3-dimensional 
orientation vector. However, this representation is not unique. One minimal representation 
that is commonly used to represent the attitude is Euler angles. The 3x1 Euler angle vector 
7] is given by [14] : £(77)77 = w, where E(r]) is a 3 x 3 transformation matrix. E(rj) becomes 
singular for certain values of tj ; however, it is to be noted that the limitations imposed on 
the allowable orientations due to this singularity are purely mathematical in nature and 
have no physical significance. The problem of singularity in 3-parameter representation 
of attitude has been studied in detail in the literature. An effective way of overcoming the 
singularity problem is to use the quaternion formulation (see [15]- [17]). 

The unit quaternion a is defined as follows. 

a = { q T , < 24 }^, a — 

a = (ai,d2,a3) r is the unit vector along the eigen-axis of rotation and <j> is the magnitude 
of rotation. The quaternion is also subjected to the norm constraint: 

a T a + ~ 1 (6) 


Qfi 

<* 3 


• A\ ,<t ) \ 

<X4 = C0S{~) 


(5) 
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It can be also shown [18] that the quaternion obeys the following kinematic differential 
equations. 

a = i(u> x a + ar4w) (7) 


1 T— 
a 4 = a 


( 8 ) 


The attitude control of a single-body rigid spacecraft using a quaternion feedback has been 
thoroughly investigated [12], [15-17]. We shall use quaternion representation for the central 
body attitude. The quaternion can be computed [18] using Euler angle measurements (Eq. 

4). 

Defining 0 = (a 4 - 1) and denoting p = z, equations (2), (7), and (8) can be rewritten as: 


Mi + Cz + Dz + K{0i x3 ,0 T ,q T ) T = B r u 


(9) 


a = x a + (0 + l)w) (10) 

= (in 

In equation (9) the matrices M and C are functions of p, and (p, p), respectively. It is to 
be noted that the first three elements of p associated with the orientation of central body 
can be fully described by the unit quaternion. Hence, M and C are implicit functions of a, 
and therefore, the system represented by equations (9)- (11) is time-invariant and can be 
expressed in the state-space form as follows: 


x - f(x,u) 


( 12 ) 


where x = (a T , /3,0 T ,q T ,z T ) T . Note that the dimension of x is (2 n + 1), which is one more than 
the dimension of system in (2). However, one constraint (Eq. 6) is now present. It can be 
easily verified from (6)-(8) that the constraint (6) is satisfied for all t > 0 if it is satisfied at 
t = 0 . 
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Stability with Dissipative Control Law 


Consider the control law u, given by: 


U — Gpf) — GrVr 


( 13 ) 


where p = {«? ,0 T } T . Matrices G p and G r are symmetric positive definite ( k x k) matrices and 
G p is given by: 


G p - 


(i + is+n )g p , o 3x(fc _ 3) 


(14) 


°(k- 3)X3 G p2(fc-3)x(fc-3) J 

Note that eqs. (13) and (14) represent a nonlinear control law. If G p and G r are symmetric 
and positive definite, this control law can be shown to render the time-rate of change of 
the system’s energy negative along all trajectories; i.e., it is a dissipative control law. The 
closed-loop equilibrium solution can be obtained by equating all the derivatives to zero in 
Eqs. (2), (10), and (11), with the input as in (13) and (14). It can be easily verified that 
the equilibrium solutions of the closed-loop system given by eqs. ( 12 ) and (13) are: a = 0 , 

0 = 0 , q = 0 , z = 0 , and 0 = 0 or -2 (i.e., a 4 = ±l). Thus, there appear to be two closed-loop 
equilibrium points corresponding to 0 = 0 (a 4 = l) and 0 = -2 (<* 4 = -l) (all the other state 
variables being zero). However, from Eq. (5), 0 - 0 (<* 4 = 1 ) => <j> = 0 , and 0 = -2 (» 4 = - 1 ) => 
<j> = 2tt, i.e., there is only one equilibrium point in the physical space. One objective of the 
control law is to transfer the state of the system from one orientation (equilibrium) position 
to another orientation. Without loss of generality, the target orientation can be defined to 
be the origin ( x = 0 ), and the initial orientation, given by (a(O),/?(O), 0 (O)) can always be 
defined in such a way that |0,(O)| < 7 r, and -1 < /?(0) < 0, i.e., 0 < a 4 (0) < 1 (corresponding 
to |^| < tt) and (a( 0 ),a 4 ( 0 )) satisfy Eq. ( 6 ). 

The following theorem establishes the global asymptotic stability of the physical 
equilibrium state of the system. 

Theorem 1. Suppose <jp 2 (^_ 3 ) x (^_ 3 ) and G r (ky.k) are symmetric and positive definite, and 
G p \ = p /3 where p > 0 . Then, the closed-loop system given by equations (12) and (13) is 
globally asymptotically stable. 
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Proof. 


Consider the Lyapunov function 

V = \p T M(p)p+ \q T kq+±6 T G p2 e+\a T (G pl +2pI 3 )a+tf 2 (15) 

V is clearly positive definite and radially unbounded with respect to the state vector 
{a T ,0,9 T ,q T ,p T } T since M(p), A, G p \ and G p2 are positive definite symmetric matrices. 
Taking the time derivative, we have: 

V = p T Mp + ip T Mp3- q T Kq + 9 T G p2 9 + a T (G p \ + 2pl 3 )a + 2p00 (16) 

Using (2), (4), (10), (11) and (14), we get: 

V = pT u p? (-M — C)p — i? Dp — p^Kp + q^Kq + 9^ G p2 9 + —{Qa)^G p \Q + — (/? 4- 1 G p \a 4- pu>^ a 

(17) 

where ft = (wx) is a skew-symmetric matrix. Substituting for u and noting that, p T I<p = 
q T kq, (Sla) T G p ia = 0, and using Property S of the system, we obtain 

V = -p t {D + B T GrB)p-{Bp) T G p p+ i(/?+ l)u T G pl a + tiuj T a + 9 T G p2 6 (18) 

Note that ( Bp) T G p p = + l)w T G pl a 4- pw T a + 9 T G p2 9. After several cancellations, we get 

V = -p T (D + B T G T B)p (19) 


Since ( D + B T G r B ) is a positive definite symmetric matrix, V < 0, i.e., V is negative 
semidefinite, and P = 0=>p = 0=>p = 0. Substituting in the closed-loop equation we get 


-B T G p p = 


-G p p 


Ofcxl 

.0(n-fc)xl . 


. kq . 


( 20 ) 


=> p = 0, and 9 = 0, i.e., a = 0,9 = 0, and (3 = 0 or -2. Thus, V < 0 along all trajectories, 
and V = 0 at the two equilibrium points. Therefore, if the system’s initial condition lies 
anywhere in the state space except at the equilibrium point corresponding to /? = -2, then 
the trajectory will asymptotically approach the origin, i.e., x = 0; and if the system is at 
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the equilibrium point corresponding to /? = -2 at t = 0 then it will stay there for all t > 0. 
However, consistent with the previous discussion, the two equilibrium points in the state 
space represent the same equilibrium point in the physical space; hence it can be concluded 
that the system is globally asymptotically stable. ■ 

A Special Case: 

Consider a special case where the central body attitude motion is small. This can occur 
in many realistic situations. For example, in the case of a space station-based or Shuttle- 
based manipulator, the inertia of the base (central body) is much larger than that of any 
manipulator link or payload. In such cases the rotational motion of the base can be assumed 
to be in the linear region, although the payloads (or links) attached to it can undergo 
large rotational and translational motions and nonlinear dynamic loading due to Coriolis 
and centripetal accelerations. For this case, the attitude of the central body is simply the 
integral of the inertial angular velocity and the use of quaternions is not necessary. The 
equations of motion (2) can now be expressed in the state-space form simply as: 


x = g(x, u) 


( 21 ) 


where x = ( p T ,p T ) T . 

The dissipative control law u is now given by: 

u - ~G p y p - G r y r (22) 

where, G p is symmetric positive definite (k x k) matrix, 

y v = Bp and y T - Bp (23) 

y p and y r are measured angular position and rate vectors. 

Theorem 2. Suppose G Pkxk and G rkxk are symmetric and positive definite. Then, the 
closed-loop system given by equations (21), (22) and (23) is globally asymptotically stable. 
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Proof. 


Consider the Lyapunov function 

V = lp T M(p)p+ \p T {K + B T G p B)p (24) 

V is clearly positive definite since M(p) and ( K + B T G V B) are positive definite symmetric 
matrices. Taking the time derivative, letting K = (K + B T G P B), and simplifying, we get 

V = p T (^M - C)p - p T Kp + i> T Kp -P T (D+ B t G t B)p (25) 

Again, using Property 5, we get, p t {\M - C)p = 0, and after some cancellations, we obtain 

V = -p T {D + B T G r B)p (26) 

Since ( D + B t G t B) is the positive definite symmetric matrix, V < 0, i.e., V is negative 
semidefinite in p and p and V = 0=>p = 0^>p = Q. Substituting in the closed-loop equation 
we get 

(K + B T G p B)p = 0 = 0 (27) 

Thus, V is not zero along any trajectories; then, by LaSalle’s theorem, the system is globally 
asymptotically stable. ■ 

The significance of the two results presented above is that any nonlinear multibody system 
belonging to these classes can be robustly stabilized with dissipative control laws. In the case of manipulators, 
this means that one can accomplish any terminal angular position from any initial position 
with guaranteed asymptotic stability. 

Robustness to Actuator/Sensor Nonlinearities 

Theorem 2 proves global asymptotic stability for the practically important case where 
the central body motion is in the linear range and the other bodies undergo nonlinear 
motion. It assumes linear actuators and sensors. In practice, however, the actuators and 
sensors have nonlinearities. The following theorem extends the results of [2] to the case 
of nonlinear flexible multibody systems. That is, the robust stability property of the 
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dissipative controller is proved to hold in the presence of a wide class of actuator/sensor 
nonlinearities. 

Let , V'pd ), and Vvt(-) denote the nonlinearities in the ith actuator, position sensor, 
and rate sensor channels, respectively. Assuming G p and G r are diagonal, the actual input is 
given by: 

«. = i>ai[-Gpill>pi(y p i) - G ri Vv«(l/rt)] (t = 1,2,..,*) (28) 

We assume that ip p i, V'a* and Vv« (i = 1,2, ...,k) are continuous single- valued functions: 

R — ► R. [A function ip(v) is said to belong to the (0,oo) sector (Figure 3) if V’(0) = 0 and 
izV’(t') > 0 for v ^ 0: rp is said to belong to the [0, oo) sector if vrp(v) > 0]. The following theorem 
gives sufficient conditions for stability. 

Theorem 3. Consider the closed-loop system given by (21), (22), (23), and (28), where 
G p and G T are diagonal with positive entries. Suppose (for * = 1,2, ...,k) y ai , V’pi, and tp r { are 
single- valued, time invariant continuous functions belonging to the (0,oo) sector and rp ai are 
monotonically nondecreasing. Under these conditions, the closed-loop system is globally 
asymptotically stable. 

Proof. 

(The proof closely follows [2].) Let w = -y p (fc-vector). Define 


't’vM) = “tW-") (29) 

i’rM = -i’ri(-i') (30) 


If ip p i, ip r i € (0,oo) or [0,oo) sector then xp pi , xp ri also belong to the same sector. Now, consider 
the following Lure-Postnikov Lyapunov function : 

j U 

i i K rw x 

V = -p T M(p)p + -q T Kq + £ / rP^Gp^u)}^ (31 ) 

.=1 Jo 

where, K is the symmetric positive definite part of K. Taking the time derivative and using 


( 2 ), 


V = p T [B T u 



(32) 
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Upon several cancellations and using Property S, 

k k 

V = Y UiVr ' - "’tV’oilGpjV’p, •(«',•)} ( 33 ) 

1=1 i=i 

where, matrix D is the positive definite part of D. 

k 

V = -q T Dq - Y '* > i('l>ai[Gri4’ri(™i) + G pi 1 p pi (Wi)] - i> ai [G pi xl> pi (wi)]) (34) 

1=1 

If tp a i are monotonic nondecreasing and ip r { belong to the (0,oo) sector, V < 0, and it can be 
concluded that the system is at least Lyapunov-stable. Now we will prove that in fact the 
system is globally asymptotically stable. First, let us consider a special case when xj) ax are 
monotonic increasing. Then V < - q T Dq , and V = 0 only when q = 0 and w = 0, which implies 
p = 0 p = 0. Substituting in the closed-loop equation, 


Kp — B? xl) a [— Gpipp(yp)\ 


(35) 


' 0 ' 


V'o {—G p ifi p (y p )} 



0 


(36) 


=> ^a[-G p xpp(y p )] = 0, and q = 0 


If \l> p i belong to the (0, oo) sector, = %i{^) = 0 only when v - 0. Therefore, y p = 0. Thus, 

V = 0 only at the origin, and the system is globally asymptotically stable. 

In the case when actuator nonlinearities are of the monotonic nondecreasing type (such 
as saturation nonlinearity), V can be 0 even if w ^ 0. Figure 4 shows a monotonically 
nondecreasing nonlinearity. However, we will show that every system trajectory along which 

V = 0, has to go to the origin asymptotically. When w ± 0, V = 0 only when all actuators 
are locally saturated. Then, from the equations of motion, it means that system trajectories 
will go unbounded which is not possible since we have already proved that the system is 
Lyapunov-stable. Hence, system trajectories have to approach the origin asymptotically, 
and the system is globally asymptotically stable. ■ 
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Concluding Remarks 

Stability of a class of nonlinear multibody flexible space systems was considered using 
a class of dissipative control laws. It was shown that robust global asymptotic stability 
can be obtained using a nonlinear feedback of the central body quaternion angles, relative 
body angles, and angular velocities. For the practically important special case wherein the 
central body motion is in the linear range, it was shown that global asymptotic stability is 
obtained with a linear dissipative control law. Furthermore, it was shown that the robust 
stability is preserved in the presense of a wide class of actuator and sensor nonlinearities. 

All the stability results presented are valid in spite of modeling errors and parametric 
uncertainties. The results have a significant practical value since the mathematical models 
of such systems usually have substantial inaccuracies, and the actuation and sensing devices 
have nonlinearities. 
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Branches 



Figure 1. Multibody system 



Figure 2. Multibody system with a single chain 


16 



Figure 3. Nonlinearity belonging to (0, oo) sector 



Figure 4. Monotonically non -decreasing nonlinearity 
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